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Abstract 

The high mass neutral quantum states envisaged by theories of physics beyond the standard 
model can at the hadron colliders reveal themselves through their decay into a pair of photons. 
Once such a peak in the diphoton invariant mass distribution is discovered, the determination 
of its spin through the distinctive photon angular distributions is needed in order to identify 
the associated nonstandard dynamics. We here discuss the discrimination of the spin-2 Randall- 
Sundrum graviton excitation against the hypothesis of a spin-0 exchange giving the same number 
of events under the peak, by means of the angular analysis applied to resonant diphoton events 
expected to be observed at the LHC. The spin-0 hypothesis is modelled by an effective interaction 
of a high mass gauge singlet scalar particle interacting with the standard model fields. The basic 
observable of our analysis is the symmetrically integrated angular asymmetry Acs, calculated for 
both graviton and scalar s-channel exchanges to next-to-leading order in QCD. 
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I. INTRODUCTION 



Diphoton final states represent a very important testing ground for the standard model 
(SM), for example they may be one of the main discovery channels for the Higgs boson 
searches at the CERN LHC. Moreover, similar to the case of dileptons, the inclusive pro- 
duction of two-photon high mass resonance states at the LHC: 

p + p->77 + ^, (1) 

is considered as a powerful, clean test of New Physics (NP), would an excess of 77 events 
be observed with respect to the prediction from the SM cross section. 

One NP scenario of particular importance is the case of the spin-2 Kaluza-Klein (KK) 
graviton excitations predicted by the Randall-Sundrum (RS) model of gravity in one warped 
spatial extra dimension This model suggests a rich phenomenology that includes the 
production of diphoton resonances, to be explored at collider energies, see, for example, 
Refs. (2 - 4 ] . The existence of such graviton excitations can be signalled by the occurrence 
of peaks in the invariant mass distribution of the photon pairs and, indeed, the lowest lying 
predicted diphoton peak has recently been searched for in experiments at the p p Fermilab 



Tevatron collider 



6[, and at the 7 TeV p p LHC collider with time- integrated luminosity 



of the order of 40 pb 1 



In these experiments, exclusion mass limits on the lightest RS 
resonance of the TeV order have been set, and graviton mass scales larger than 1 TeV will 
certainly be in the kinematical reach of LHC. 

Assuming that a diphoton peak at an invariant mass value Mr is observed, its association 
to a specific NP scenario would be possible only if we are able to discard other competitor 
models, potential sources of the peak itself with same Mr and same number of events. 
Basically, for any nonstandard model one can define, on the basis of the foreseeable statistics 
and uncertainties, a discovery reach on the relevant heavy resonance R as the upper limit 
of the range in Mr where, in a specific domain of the model parameters called "signature 
space" , the peak is expected to give a signal observable over the SM prediction to a prescribed 
confidence level. Instead, the identification reach on the model is the upper limit of the range 
in Mr where it can be identified as the source of the peak, once discovered, or, equivalently, 
the other competitor models can be excluded for all values of their respective parameters. 
Of course, for many models, identification should be possible only in a subdomain of their 
signature space. 



The determination of the spin of an observed resonance clearly represents an important 
selection among different classes of nonstandard interactions. In the case of the inclusive 
diphoton production ([I]), the tool to directly test the spin-2 of the RS graviton resonance or, 
equivalently, exclude the hypothesis of a spin-0 scalar particle exchange, would be provided 
by the distinctive angular distributions in the angle 9 between the incident quark or gluon 
and the final photon in the diphoton center-of-mass frame. This is similar to dilepton 
production, the difference being that in this case the hypotheses of both the spin-0 and the 
spin-1 exchanges must simultaneously be excluded. 

The spin-2 test of the lowest-lying RS graviton in lepton-pair collider events, through the 
direct comparison of the angular distributions for the various spin hypotheses, was earlier 
discussed in several papers, see, e .g., Refs. j9-ll|, and experimental angular analysis were 



attempted at the Tevatron in Ref. 



12| . A potential difficulty of the direct- fit angular analysis 



at the LHC is that generally, due to the symmetry of the proton-proton initial configuration, 
the determination on an event-by-event basis of the direction of the initial parton, hence of 
the sign of cos 9, is in principle not fully unambiguous, so that cuts in phase space must be 
applied in this regard. 



The spin-2 RS graviton analysis of LHC dilepton events proposed in Ref. 



131 ] . makes use 



of a "center-edge" angular asymmetry Ace where the above mentioned ambiguity should 



not be present fl^ . [15|. Essentially, in this observable the dilepton events are weighted 
according to the cos 9 differential distributions, and the asymmetry is defined between cross 
sections symmetrically integrated over "center" and "edge" angular intervals. Recently, 
asymmetries conceptually analogous to Ace, have been applied to heavy quantum states spin 
identification in Refs. 



16 



171 ]. and a comparison of the performances of different methods 



for heavy resonances identification has been presented in Ref. [18]. Angular analyses for 
different spin-mediated Drell-Yan processes have been applied to a variety of NP models in 
Ref. [19]. 

Here, we propose the application of Ace to the angular analysis of the diphoton produc- 
tion process ([I]) at the LHC. As remarked previously, the selection of the spin-2 RS graviton 
amounts in practice to exclude the hypothesis of a spin-0 particle exchange with same mass 
Mr and producing the same number of diphoton events. Ideally, one advantage of the 
diphoton channel over dileptons can be represented by the doubled statistics expected in 



the former case 



201 ] . Also, the automatic exclusion of the spin-1 hypothesis [2l|, |22j, should 
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in any case allow a simplification of the analysis from the phenomenological point of view. 
Finally, the consideration of process ([1]), in addition to dilepton production, is needed for 
an exhaustive test of model [1]. 

For our analysis we have used the calculations of the required differential cross sections 
to next-to-leading order (NLO) in QCD, and this is essential at a hadron collider as the 
theoretical uncertainties get reduced when higher order corrections are included. Further- 
more, as a result of new interactions in a NP model, there will be additional subprocesses 
that contribute at leading order (LO) itself (e.g., gg — > 77 in the RS model) and hence the 
signal can receive enhanced contributions due to the NLO corrections. 

Specifically, in Sec. II we review the definitions of the basic cross sections involved in the 
asymmetry Aqe] Sec. Ill will be devoted to the relevant properties and the characteristic 
angular distributions for the RS graviton and for the competitor scalar particle exchanges 
in process (pQ), for which we will adopt the model recently proposed in Ref. {23]. In Sec. 
IV we discuss the NLO QCD effects to the diphoton production rates and to the angular 
distributions, for both kinds of spin exchange. Sec. V contains an outline of the AcE-based 
angular analysis and the consequent numerical results for RS identification, in the LHC 
center-of-mass running configurations y/s = 14 TeV and \/s = 7 TeV. Finally, Sec. VI 
contains some conclusive remarks. 



II. CROSS SECTIONS AND CENTER-EDGE ASYMMETRY 

The total cross section for a heavy resonance discovery in the events (pQ) at a diphoton 
invariant mass M = Mr can be expressed as 

^ ll) = J_J z J M ^ iM —, (2) 

where the rapidity of the individual photon |?7 7 | < 2.5 and z = cos 9 is chosen such that 
\z\ < 0.98. 

Resonance spin-diagnosis uses the comparison between the characteristic photon differen- 
tial distributions for the two hypotheses for the resonance bump, R = G spin-2 Kaluza-Klein 
(KK) modes of the RS model graviton and R = S a massive scalar: 

da f Mn+AM/2 da 

dM——. (3) 



dz J Mr -am/2 dMdz 



In Eqs. ([2]) and (J3J), cuts on the phase space accounting for detector acceptance are implicit, 
and AM is an invariant mass bin around Mr, which should somehow reflect the detec- 
tor energy resolution and be sufficiently large as to include the resonance width. In the 



calculations worked out in the sequel, we use for this mass window the expression 

AM = 24 (0.625M + M 2 + 0.0056) 1/2 GeV. (4) 

Actually, Eq.fll]) was derived in connection to the ATLAS and CMS experiments on dilepton 
production, but we assume it also for the calculations of diphoton production of interest here. 
Obviously, for a resonance sufficiently narrow, the integral over M should be practically 
insensitive to the size of AM, whereas it should be essentially proportional to AM for a 
flat background such as the SM. Besides |r/ 7 | < 2.5 mentioned above, the assumed typical 
cuts on harder (softer) photons are p]_ > 40(25) GeV, and the statistics will be estimated 
by taking a photon reconstruction efficiency e 7 = 0.80. 

Moreover, to evaluate Eqs. (j2J) and (j3J), the partonic cross sections will be convoluted 
with the CTEQ6L and CTEQ6M parton distributions sets for LO and NLO cross sections, 



respectively, with Aqcd — 0.226 GeV [25j. In particular, for resonance discovery, process ([I]) 
must be observed with a number of events well-above the background from SM processes. 
Specifically, denoting by N$ and iVgM the numbers of signal and SM events in the 77 invariant 
mass window, the statistical significance of a 5-cx signal would be ensured by the criterion 



that Ns should be larger than max(5-y/iVsM, 10). 

The z-evenly-integrated center-edge angular asymmetry A C e is defined as: 

da 



A C e = — with a CE 
a 



dz. (5) 

dz 



In Eq. (J5J, < z* < z cut defines the separation between the "center" (\z\ < z*) and the 
"edge" (z* < \z\ < z C ut) angular regions and is a priori arbitrary to some extent. In previous 



applications, see for example Refs. 



13 



3Hl7j. the "optimal" numerical value turned out to be 
z* ~ 0.5, and we shall keep this value of z* here as well. One can notice that by definition 
Aqe is symmetric under z -H- —z, hence it is insensitive to the sign of z. Moreover, as being 
a ratio of integrated cross sections, an advantage of Aqe is that it should be less sensitive to 
theoretical systematic uncertainties, such as the uncertainties from different sets of parton 
distributions and from the particular choice of factorization and renormalization scales. 
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III. GRAVITON RESONANCE AND SCALAR EXCHANGES 



We here sketch the models we are interested in, together with their features relevant to 
the resonance spin and the distinctive angular distributions for process (prj). 

A. RS model of gravity with one compactified extra dimension 

This model, originally proposed as a solution to the gauge hierarchy problem Mew ^ 
Mpi, consists of two 3-branes, and one compactified warped extra spatial dimension y with 
exponential warp factor exp (— kn\y\) [1]. Here, k > is the 5D curvature, assumed to 
be of the order of Mpi. The two branes are placed at orbifold fixed points, y = with 
positive tension called the Planck brane and the second brane at y — R c with negative 
tension called the TeV brane. The basic, simplifying, hypothesis is that the SM fields are 
localized on the TeV brane, whereas gravity originates on the Planck brane but is allowed 
to propagate everywhere in the 5D space. The consequence of this setup is the existence of 
KK modes of the graviton that can be exchanged in the interactions among SM particles 
in TeV brane. Owing to the exponentially suppressing warp factor, mass scales, in passing 
from the Planck brane to the TeV brane, can get the size of the TeV order. Moreover, a 
specific mass spectrum of such KK resonances is predicted in terms of an effective mass 
scale defined as = Mpjexp (—kirR c ), that for kR c ~ 12 happens to be of the TeV order 



(here, Mp\ = 1/v8vtGn with the Newton constant). These resonances, represented by 
spin-2 fields hfy \ can in process ([I]) show up as (narrow) peaks in M 77 = M, through the 
interaction 

c= T »vJ2h$. (6) 

iW n=l 

Here: T^ v is the energy- momentum tensor of the SM; the characteristic mass spectrum is 
M n = x n kexp (— kiiR c ) (with x n the roots of the Bessel function J\{x n ) = 0); and the 
resonance widths are T n = pM n x^(k/ Mpi) 2 , with p a calculable constant depending on the 
number of open decay channels, of the order of 0.1. 

The model can therefore conveniently be parametrized in terms of Mq = Mi, the mass 
of the lowest graviton excitation, and of the "universal" dimensionless graviton coupling 
c = k/Mp\. Theoretically, the expected "natural" ranges for these parameters, avoiding 
additional mass hierarchies, are: 0.01 < c < 0.1 and < 10 TeV |2|. The 95% CL 



6 



experimental lower bounds on Mq from previous analysis vary, essentially, from 0.6 to 1.0 
TeV as c ranges from 0.01 to 0.1 , l26j. Quite recently, preliminary results from RS 
graviton searches in dilepton inclusive production at the 7 TeV LHC with luminosity 1.2 
fir 1 , indjcate,95% CL lower limits on M G of 0.7 TeV for c = 0.01 up to 1.6-1.7 TeV for 
c = 0.1 



27 



28|. 



In hadronic collisions, in QCD at LO, photon pairs can be produced via the quark- 
antiquark annihilation q + q — > 7 + 7, and the gluon-gluon fusion g + g — > 7 + 7. The 
relevant diagrams at this order, for the SM and the RS graviton exchange, are represented 
in Figs. [TJ and [2J Actually, the SM box diagram in Fig. [2] is of higher order in a s and, as 
discussed in the next section, for the values of the 77 invariant mass M in the TeV ran ge o f 



interest here, its contribution turns out to be negligible, as earlier noticed also in Ref. 



29). 



(WW 7 



Vwv 7 





FIG. 1: Feynman diagrams contributing to the subprocess qq —> 77, where the dashed line rep- 
resents either a RS graviton KK mode G or a scalar exchange S. The crossed diagrams are not 
displayed. 




FIG. 2: Feynman diagrams contributing to the subprocess gg —> 77, including the exchange of 
either a RS graviton KK mode G or a scalar S which is denoted by the dashed line. The crossed 
diagrams are not displayed. 

Using Feynman rules for graviton exchange {20I, [jl]], the z-even angular dependencies 



7 



needed in (jSj) can be written as 3lM33j: 
da(qq — > 77) 1 



dz 



192 vrs 

- 47T«Q2 S 2 R e [C(s)] {l + Z 2 ) 



+ ^|c(s)l 2 (i-/) 



(7) 



where the subprocess Mandelstam variable is the diphoton invariant mass, v£ = M, and 
a is the electromagnetic coupling constant with Q q the quark electric charge. Moreover, in 
Eq. (J7J), C(s) represents the sum of KK graviton propagators with masses M = M n and 
widths r„: 



(8) 



Al ^ s-M* + iM n T n 

n " 

In practice, from the phenomenology, just the lowest graviton mass Mq = Mi and, perhaps, 
the next one M2 at most, can be expected to fall within the discovery reach of LHC. 

The cross section for the gg —> 77 subprocess via the RS graviton excitation exchange is 



d(r{gg -»■ 77) 



|C(s)| 2 (l + 6^ + ^ 4 ). 



(9) 



dz 8192 Ti 

Notice that a factor 1/2 is embodied in Eqs. (J7|) and (Q to account for the identical final 
state photons. 



B. The model for scalar exchange 

In principle, in order to discriminate the graviton spin-2 angular distribution from the 
spin-0 hypothesis, we might limit ourselves to make a comparison of Eqs. (JTj) and ([9]) with the 
results of a generic flat (in z) distribution numerically tuned to the same number of events 
around Mq. However, for the graviton exchange we shall use a description supplemented 
by the cross sections calculated to NLO in QCD and, to consistently fully exploit the QCD 
dynamics also for the spin-0 scenario, we need an explicit model for a scalar particle exchange 
with definite couplings to the SM fields, in particular to photon pairs. 



We consider the simple model of a scalar particle S, sing 
and with mass M = Ms of the TeV order, proposed in Ref. 



et under the SM gauge group 



23j. The trilinear couplings of 



S with gluons, electroweak gauge bosons and fermions, are in this model: 

£ = c, 9 fG%G a + c 2 9 ^w; y w i + Cl 9 -1 B , U B^S + Cf^ffS. (10) 
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In Eq. ( ITOj) . A is a high mass scale, of the TeV order of magnitude, and c's are dimensionless 
coefficients that are assumed to be of order unity, reminiscent of a strong novel interaction. 
In our subsequent analysis we take, following Ref. 23|, A = 3 TeV. As for the c's, we shall 
leave their numerical values free to vary in a range of the order of (or less than) unity, but 
constrained so that the scalar particle width Tg could be comparable to (or included in) the 
diphoton mass window AM of Eq. fl4]). 

The leading order diphoton production process is in this model dominated by the s- 
channel exchange gg — >■ S — >■ 77. As it is of order a 2 at LO, it will be sensitive to the choice 
of this coupling constant. The Feynman diagrams in the scalar model will be similar to 
those in the RS model, as shown in Figs. Q] and Fig. [2] with the KK mode replaced by the 
scalar. The corresponding differential cross section reads, including a factor 1/2 for identical 
final particles: 

da{gg 77) _ 1 1 f c 3 g 2 \ 2 ( ( Cl + c 2 )e 2 \ 2 §3 \ D ^2 ^ 



dz 216tt V A / V A 

In Eq. ([HJ, D(s) is the scalar propagator, 



D{§) = l-W + iMsTs ' <12) 



and the expression of the total width Tg in terms of the c's and A introduced in Eq. ( flOl) 
can be obtained from 23], by summing the partial widths reported there into gg (dominant, 
as is the cross section, proportional to c 3 ), 77, W + W~, ZZ, and //. 

In Fig. [3] we represent, as an example, the domains in C3 vs all other parameters q (i ^ 3) 
assumed equal to each other, allowed by the constraint Tg < AM for different values of Ms- 



As the coupling of the scalar particle to quarks is proportional to the quark mass much 
smaller than A, one can neglect at LO the subprocess qq — > S — > 77, and the cross section 
for pp — > S — > 77 will be expressed simply by the gluon-initiated process plus the SM 
contribution appearing in Eq. Thus, at this LO, there is in the considered spin-0 

model no interference between the SM quark-initiated contribution and the scalar-exchange 
amplitude. However, quark contributions cannot be neglected at the NLO QCD order and, 
in particular, a (small) interference of the gg SM box diagram with the gg scalar exchange 
will occur, as it will be specified in the next section. 

Preliminary to the numerical analysis of A^, in the next section we briefly describe the 
estimate of the next-to-leading order QCD effects. 

9 



I I I I I I I I I I I I I I I I I I I I I I I I I- 

|- M s = 1 TeV 



0.7 
0.5 
0.3 
0.1 



M s = 1.5 TeV 
M s = 3 TeV 



- i . . i . . \ ■ i ■ ■ i . . \ 



I . . I . . 'U 



0.1 0.4 0.7 1. 1.3 1.6 1.9 2.2 2.5 



FIG. 3: The dashed lines in the parameter plane (cj,C3) with equal Ci (i ^ 3), correspond to the 
equation Ts = AM for Ms = 1 TeV, 1.5 TeV and 3 TeV, A = 3 TeV. The allowed area is at the 
left of these lines. 

IV. NEXT-TO-LEADING ORDER QCD EFFECTS 

The diphotons produced in a hadronic collision could originate from the hard partonic 
interaction (direct photon) or at least one photon could be produced in the hadronisation 
of a parton (fragmentation photon). At higher orders in QCD, there would be final state 
collinear singularities of QED origin as a result of the emission of a photon from a quark. 
These singularities can be factored out and absorbed into the fragmentation functions. The 
fragmentation functions are additional non perturbative inputs that are not well-understooc . 
An alternate approach to isolate direct photons is the smooth cone isolation criterion 34J, 
which ensures that the fragmentation contributions are suppressed without affecting the 
cancellation of the conventional QCD singularities. 

We start from the extra dimension scenario, where the NLO corrections in QCD were con- 



sidered for the phenomenologically interesting process like the dilepton 
ton production 



35 



38 



- 1371 ] and dipho- 



at hadron colliders. The essential Feynman rules for KK modes 



coupling to ghosts and KK mode propagator in (4 + e)-dimensions needed for the NLO 
computation have been introduced in [35 . For the diphoton computation a semi-analytical 
two-cutoff phase space slicing method [40( to deal with various singularities of infrared (IR) 
and collinear origin that appear at NLO in QCD was used. After cancellation and mass 
factorisation of these singularities, the remaining finite part is numerically integrated over 
the phase space by using Monte-Carlo techniques. 

As stated in previous sections, diphotons can in the SM be produced to LO in QCD, 
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in the quark anti-quark annihilation subprocess qq — > 77. Photon pairs produced via the 
gluon fusion subprocess through a quark box diagram gg — > 77, though of the order a^, 
have cross sections comparable to those of the qq — > 77 subprocess for the low diphoton 
invariant mass. In the light Higgs boson searches, this subprocess plays an important role, 
due to the large gluon flux at small fractional momentum x, and is formally treated as a 
LO contribution although it is of O(a^), hence is in reality a next-to-next-to leading order 
contribution. However, it falls off rapidly with increasing diphoton invariant mass and in 
the mass range of interest for the TeV scale gravity models it need not be included at LO. 



It has been demonstrated in 



38 



that the contribution of this subprocess in the SM is 
few orders of magnitude smaller than that of the qq sub process for diphoton invariant mass 
M > 500 GeV. 1 



SCOOT? 
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TOOT?5 





FIG. 4: Feynman diagrams contributing to the subprocess qq — > 775. The dashed line could be 
either R = G or S depending on the model. For the scalar exchange S the last diagram involving 
a four point vertex does not exist. 



The NLO QCD corrections to the diphoton production process would involve real emission 
diagrams corresponding to the following 2 — )■ 3 subprocesses: (a) qq —> 77^ (Fig. H]), (b) 
gg — > 77# (Fig. [5]) and (c) gq(q) —> 77Q'(g) (the gq diagrams can be obtained from the qq 



1 See also the discussion of the SM NLO predictions for process (TT]) in Rcf. 41 1. 
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1 




FIG. 5: Feynman diagrams contributing to the subprocess gg — > 775, where the dashed line 
corresponds to a spin-2 RS KK mode G or a scalar exchange S. In the case of a scalar exchange 
the last diagram involving a four point vertex does not exist. 



diagram by appropriate inter change of initial and final state particles). In the virtual part, 
0(a. s ) corrections at one loop arise as a result of the interference between the one loop graphs 
at 0(a 8 ) of (SM + NP) and the Born graphs at O(a s ) of (SM + NP). Some of the typical 
NP one loop Feynman diagrams to 0(a s ) are shown in Fig. |6j The qq channel gets such 
contributions from both the SM and the graviton exchange, while in the gg channel the SM 
contribution already begins at 0(a s ). The SM gg — > 77 subprocess amplitude can interfere 
with the gluon initiated LO gravity exchange subprocess, giving a 0(a s ) contribution Fig. 
|2J In the qg channel there is no virtual contribution to this order in either the SM or the 
NP models of interest here. 

Both the virtual and real corrections have been evaluated with 5 quark flavors and in the 
limit of vanishing quark masses. The n-point tensor integrals appearing from integration over 
loop momenta were simplified using the Passarino-Veltman reduction 42j, the computation 



was done using dimensional regularization in (4 + e)-dimensions, in the MS scheme. 
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FIG. 6: Feynman diagrams contributing to one loop subprocess to the qq and gg channel in the 
NP models. 

The virtual contribution here does not contain UV singularities, this can be attributed 
to the facts that i) the electromagnetic current is conserved and does not receive any QCD 
corrections, and ii) that the graviton couples to the energy-momentum tensor of the SM 
fields which is a conserved quantity and therefore is not renormalized. The poles in e arise 
from loop integrals and correspond to the soft and collinear divergences, configurations where 
a virtual gluon momentum goes to zero give soft singularities while the collinear singularities 
arise when two massless partons become collinear to each other. 

The 3-body phase space of the real emission diagrams has regions which are soft and 
collinear divergent. The phase space can accordingly be separated into soft (s) and hard 
regions and, furthermore, the hard region can be separated into hard collinear {HC) and 
hard non-collinear (HC) parts as follows: 

da real = da r s ea \5 s ) + dcT r $j{5 a , S c ) + da r ^(S s , 5 C ) . (13) 

The small cut-off parameters 5 S and 5 C set arbitrary boundaries for the soft (gluon energy 
< E g < 5 s ^/I/2) and collinear (0 < ty < 5 c s) regions, respectively. Here, for the 2—^3 
process with momenta pi + p 2 = P3 + Pi + P5, the Mandelstam variables are defined as 

13 



Sij = (pi +Pj) 2 , tij = (pi —pj) 2 , with s = §i2. In these mutually exclusive soft and hard 
collinear regions, the 3-body cross section can be factored into the Born, 2 — > 2, cross section 
and the remaining phase space integral can easily be evaluated in (4 + e)-dimensions to get 
the expansion of the soft and collinear singularities in powers of e. All positive powers of the 
small cut-off parameters 5 S and 5 C are set to zero, only logarithms of the cut-off parameters 
are retained. Adding the virtual and real contributions, all double and single poles of IR 
origin are automatically cancelled between the virtual and the first two terms of Eq. ( fl3l) . 
The remaining collinear poles are absorbed into the parton distribution functions. 

The hard non-collinear part in Eq. f[T3"j) corresponds to the 3-body phase space, which 
by construction is finite, and can be evaluated in 4-dimensions (e = 0). The sum of real 
and virtual contributions is now free of singularities and can hence be evaluated numerically 
using a Monte-Carlo method. It can be further seen that the explicit logarithmic dependence 
on S s and 8 C in the 2-body phase space {da r s eal (5 s ) + da r I fc(5 s , 5 C )) is cancelled by the implicit 
dependence of the 3-body hard non-collinear part da^^(S s ,5 c ) on these parameters, after 
the numerical integration is carried out. The sum of the 2-body and 3-body parts in Eq. 
(fl~3]) would be independent of the slicing parameters S s and 5 C , which is explicitly verified 



38, 



39| can be used to study the 



before the code is used for the analysis. Now these codes 
full quantitative implication of the NLO QCD corrections to the various distributions of 
interest in the extra dimension searches. 

Turning to the scalar-exchange model, the matrix elements corresponding to the inter- 
ference between the SM box diagram and the LO gg initiated tree level diagram is given 
by 



IsM*s( z ) 



9l 



^Q 2 q Re [D(s)\ 



In 



1 




Lio 



[A^ 2 - 1] 

-t 

s 

—u 



C3 9 s 

A 



(ci + c 2 ) e s 
A 



In 



-u 



t 



Li 2 



a 



2C(2) 



(14) 



where s, t = — 1(1 — z), u — — 1(1 + z) are the usual Mandelstam variables, A^ is the colour 
degrees of freedom and £(2) = ir 2 /6 is the Riemann zeta function. 

In the scalar model, the UV divergences in the virtual corrections to the gg initiated 
diphoton production process are removed by renormalization. The remaining finite contri- 
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but ion is given by 

where da/dz is the Born contribution given in Eq. ( fTTT) . The remaining soft and collinear 
finite contributions coming from the real corrections to the LO scalar model diagram are 
proportional to the Born cross section. As they originate from pure QCD, they are indepen- 
dent of the hard scattering process, hence they will be the same as those for the s-channel 
diphoton production process in the RS model and are given in Ref. 
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At the NLO in QCD, the following three subprocesses contribute to the scalar model 
cross section: i) qq — > gS ii) qg —> qS and Hi) gg — > gS, all followed by S — > 77. The 
Feynman diagrams for the qq and gg channels are given Figs. S] and respectively, wherein 
the dashed line now represents the scalar S and the four point diagrams will be absent for 
the scalar case. The couplings of S to the light quarks are proportional to the masses and 
hence are negligible so that the main contribution to the qq channel is from the S coupling 
to the gluons. Again the qg channel is related to the qq channel by suitable change of initial 
and final state. 

Fig. [7] shows, as an example, for the diphoton resonance mass values Mr = 1 and 3 TeV, 
the angular distributions at the 14 TeV LHC, at LO and NLO, in the cases of the SM, the 
RS graviton exchange with coupling constant c = 0.01, and the scalar particle exchange 
with couplings (c 3 = q = 0.3) and (c 3 = 0.3, q = 1.0), i ^ 3. 

The K-factors represent quantitatively the magnitude of the NLO QCD corrections, and 
are defined as the ratio of the NLO cross section to the corresponding LO one as follows: 

_-NLO I _NLO 

K - a LO + (J LO • ( 16 ) 
°SM ^ °NP 

Here, the subscript NP refers to the New Physics contribution (extra dimension or scalar 
exchange) and its interference with the SM to LO or NLO as the case may be, while the 
subscript SM refers to only the SM contribution. The ratio could be of total cross sections 
as is the case of Fig. [HJ or of angular distributions as is the case of Figs. [9] and [IOJ Including 
higher order QCD corrections to an observable at the hadron collider reduces the dependence 
on the factorization and renormalization scales, further since the K-factor itself could be 
large, it is essential to include these corrections. 

The dependence of the K-factors on the coupling constant c in the RS case, exhibited in 
the lower panel of Fig. [8j can be understood by the fact that for low values of c ~ 0.01 the 
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FIG. 7: Angular distributions of photons for process p + p — > r yy + X at LO and NLO QCD in SM, 
RS model with graviton coupling constant c = 0.01, and scalar resonance model with C3 = 0.3, 
Ci = 0.3, i 7^ 3 (left and right top panels) and C3 = 0.3, q = 1.0 (left and right bottom panels). 
M R = 1 TeV and 3 TeV at the 14 TeV LHC. 

NP interference with the SM could be of the same order as the pure NP part. For values 
of c > 0.05, the SM contribution is negligibly smaller than the pure NP part and can be 
neglected both in the numerator and in the denominator of Eq. ( TT6|) . which in this case is 
determined by the NP solely so that the dependence from c cancels. It is interesting to 
notice from the upper panel of Fig. |8] that the K-factors in the scalar case are larger than 
those in the RS case: the NLO corrections have enhanced the cross sections but, as shown 
in Figs. |9]and[T0l they did not noticeably change the shape of the angular distributions from 
the flat behavior of the pure scalar particle exchange cross section in Eq. ( TTTj) . 

Also, one may remark that in the example discussed here the RS graviton couples to 



quarks and gluons with equal coupling strength whereas, in the model of Ref. 23| we have 
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FIG. 8: K-factor vs M R for 14 TeV (left panel) and 7 TeV (right panel) LHC. 




FIG. 9: K-factor vs cosO for the 7 TeV LHC. 



adopted, the scalar particle couples mainly to gluons (couplings to quarks are identically zero 
in the limit of vanishing quark masses). Furthermore, the production of a scalar particle in 
the gluon-gluon fusion subprocess is qualitatively equivalent to the Higgs boson production 
in the limit of infinite top quark mass: it is well-known that K-factors (due to NLO QCD 
corrections) for the Higgs production process at hadron colliders are very high and can easily 
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FIG. 10: K-factor vs cosd for the 14 TeV LHC. 



be greater than 2.0 in the light Higgs mass region |43H45j . Hence, a similar pattern of K- 
factors can be expected in the case of a scalar production (followed by decay to photons) in 
the model considered here, see Figs. [HI El and [TUJ 
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V. ANGULAR ANALYSIS AND RS GRAVITON IDENTIFICATION 



The ylcE-based angular analysis will essentially proceed as follows. The first step will 
be the determination of graviton-scalar "confusion regions", namely, of the subdomains in 
respective discovery signature spaces of coupling constants and masses where, for Mr = 
Mq = Ms, the two models predict equal numbers of resonance signal events, N$, hence are 
not directly distinguishable on a statistical basis. In such confusion regions, one then can 
try to discriminate the models from one another by means of the different values of the Aqe 
asymmetry generated by the respective photon angular distributions. The representation of 
Aqe vs z* predicted at NLO for the two models with diphoton resonance masses of 2 and 4 
TeV, at the 14 TeV LHC, is shown in Fig. [TTJ 

NLO; c=0.01; ci=0.3; M R = 2 TeV NLO; c=0.1; ci=1.0; M R = 2 TeV 




z* z* 



FIG. 11: NLO Aqe asymmetry vs z* for the inclusive process ([T]) at the 14 TeV LHC for the 
RS model with graviton coupling constant c = 0.01 (0.1), and the scalar resonance model with 
C3 = 0.3, Cj = 0.3 (C3 = 0.3, Cj = 1.0). Mr = 2 TeV: top-left panel (top-right panel). Bottom 
panels: same as for top panels, but with Mr = 4 TeV. 
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One can start from the assumption that an observed peak at M = Mr is due to the 
lightest spin-2 RS graviton (thus, that Mr = Mq), and define a "distance" from the scalar- 
exchange model hypothesis as: 

A^ce = Aq E — Ac E . (17) 

An indication of the domain in the (Mq, c) RS parameter plane, where the competitor spin-0 
hypothesis giving same number of resonant events for Mr = Ms = Mq can be excluded by 
the starting RS graviton hypothesis, can be obtained from a simple-minded x 2 -like numerical 
procedure, similar to that used in Ref. 13[. The comparison of the deviations ( 1T7|) to the 
statistical uncertainty SAq E pertinent to the RS model, suggests the following criterion for 
spin-0 exclusion: 

X 2 = |AAg E /Mg E | 2 > X l L . (18) 

Eq. (I18p shows the definition of the \ 2 i an d Xcl specifies a desired scalar-exchange exclusion 
confidence level (for example, Xcl = 3-84 for 95% CL). With AAg E calculated in terms of 
Mr and model coupling constants, this condition will define the domain in the confusion- 
regions of model parameters where the RS spin-2 hypothesis can be discriminated from the 
scalar exchange. With (^ce) 2 m uch smaller than unity for values of z* around 0.5, we have 
to a good approximation: 



<^e^/^^,/— (19) 

V 5,min V iv 5,min 



2 

Ns min " V Nc 



where Ns, m i n will be the minimum number of RS resonance events needed to satisfy the 
criterion f[T8"j) . hence to exclude the spin-0 exchange model with same Mr = Mq = Ms in 
the confusion-region of the parameters. The knowledge of Ns, m i n determines, in turn, the 
RS resonance identification subdomain in the (Mq, c) parameter plane. 

We apply this procedure to the design LHC running conditions y/s = 14 TeV and time- 
integrated luminosity C mt = 100 fb~ . Figure H2l shows the RS graviton signature domain 
for 0.01 < c < 0.1 at LO (left panel) and to NLO (right panel) in QCD. The red line labelled 
as "Discovery" indicates the minimum number of events statistically needed for the RS KK 
resonance discovery in process at the 5-cr level, as anticipated in Sec. II. The domain 
between the lines labelled as c = 0.01 and c = 0.1 represents the number of events for RS 
KK graviton production followed by decay into a photon pair, theoretically evaluated as 
described in the previous section, for different values of Mq. The scalar resonance signature 
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FIG. 12: Number of resonance (signal) events N s vs M R (R = G, S) at the LHC with y/s = 14 TeV 
and £ m t = 100 fb -1 for the process pp — > G — > 77 + X at LO (left panel) and NLO (right panel) 
QCD. The yellow area corresponds to the KK graviton signature space for 0.01 < c < 0.1; the 
signature space for the scalar S resonance is represented by the two specific cases, C3 = q = 0.3 
and C3 = 0.3, q = 1.0. 
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space is also included in this figure, at the same LO and NLO, for simplicity by the repre- 
sentative lines labelled as C3 = q = 0.3 and C3 = 0.3, q = 1.0 (i 7^ 3), respectively. One can 
see that, for these values of the scalar coupling constants, there is a finite confusion region 
where, at given mass Mr = Mq = Ms, the numbers of predicted resonance signal events can 
be equal. Actually, such a confusion region might easily be extended to almost completely 
overlap with the full RS signature space by partially weakening the condition T$ < AM. 
This condition is, anyway, to be understood in a qualitative sense, so that more numerical 
freedom to the scalar coupling constants of Eq. (|T0|) might be allowed. Indeed, if the width 

turned out to be larger than AM in Eq. 01]), the analysis proposed here should still be 
viable, and could in this case discriminate a narrow KK graviton vs the scalar resonances, 
both by the angular observable Aqe and by the size of the widths themselves. Also, we can 
remark that, as relying just on specific angular distributions, the kind of analysis proposed 
here should be applicable more generally, to the identification of the RS graviton excitation 
from different scalar exchange models than studied here. 

The line labelled as "ID" in Fig. [12] essentially represents the solution of Eq. (1181) relating 
Mq to c, i.e., it is the minimum number of events needed for a discovered RS graviton reso- 
nance to be identified at 95% CL against the scalar particle exchange hypothesis, according 
to the AcE-based angular analysis. The differences between the left and right panels clearly 
show the need to account for the large NLO QCD effects in the theoretical description of 
the resonant diphoton inclusive production at LHC. 

In practice, as one can read from the right, next-to-leading order panel in Fig. [12], if an RS 
graviton excitation will be discovered in diphoton events at the 14 TeV LHC with 100 fb _1 
luminosity, its spin may be identified for Mq up to, roughly, 2-3 TeV for c = 0.01 — 0.1 if 
the number of observed RS resonance signal events will be larger or equal to those indicated 
by the line "ID": namely, Ns y min ~ 50 — 120. Of course, these indications follow from the 
criterion ([18]) outlined above, hence rely on statistical arguments and theoretical calculations 
of the relevant cross sections, in particular one may notice in this regard that the SM 
background turns out to be completely negligible with respect to the signal. The detailed 
assessment of the "experimental" backgrounds to the resonance discovery in process ([1]), and 
of the related systematic uncertainties, is out the scope of this paper, our purpose here is to 
just compare (and discriminate) two different theoretical explanations for same resonance 
mass events, once observed, on the basis of NLO calculations in QCD. 
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The grey area in the right panel of Fig. [12] represents the 95% CL exclusion, where the 
RS resonance should not be observable, if we account for the lower limits on Mq vs c derived 
from the dilepton production analyses recently presented in Refs. (2^] and 28]. Thus, the 
range in Mr of interest would start, in view of the new LHC results, from Mq = 1-6 — 1.7 
TeV for c = 0.1 and M G = 0.7 TeV for c = 0.01. 

Also drawn in Fig. [12] is the line A„- = 10 TeV: the "theoretical" condition A n < 10 TeV 
mentioned in Sec. III., if enforced literally, would dramatically constrain the RS discovery 
domain in the plane (Mq, c) to the events located above this line. However, also this 
condition should be understood in a qualitative sense, as is the case, in principle, of the 
assumed range of values for c. 

To make a contact with the current LHC running conditions, we repeat the RS iden- 
tification procedure outlined above for the 7 TeV case, with time-integrated luminosities 
£i nt = 1 and 10 fb -1 . The analogue of Fig. [12], but this time for the NLO calculations only, 
is represented in Fig. [THl with the same significance of the symbols. The interpretation of 
the left and right panels in this figure is also completely analogous. For example, if a RS 
resonance were discovered in diphoton events at 1 fb" 1 , its spin-2 might be identified, to 95% 
CL, up to Mq = 1.2 TeV for c = 0.1, provided the collected signal were about 60 events; 
and for 10 fb _1 , identification would be possible up to Mq = 1.5 TeV for c = 0.1 with a 
collected signal of, say, about 100 events. 

However, the situation is drastically modified by the 95% CL exclusion limits from the 



dilepton analysis of Refs. 27], [28[ at 1.2 fb" 1 , reported in Sec. Ill and represented by the 
grey areas in both panels of Fig. [13] The experimental limit Mq > 1.6 — 1.7 TeV for c = 0.1 
is not quite inconsistent with the left panel of Fig. [13] which shows that for these values the 
theoretically predicted statistics for RS events falls below that needed for 5-a discovery. The 
exclusion range starts from the "low" values c = 0.01 and Mq = 0.7 TeV, so that the left 
panel of this figure shows that, at the present stage, in principle there may be a little corner 
left available for discovery, roughly Mq between 0.7 and 1.3 TeV and small c < 0.05. On the 
other hand, this panel clearly indicates that there is no room for RS graviton identification 
at 1 fb" 1 luminosity, at least with the angular analysis presented here. Moving to the 10 
fb _1 case, the right panel of Fig. [13] shows that, with the current LHC limits, discovery 
might still be possible up to about Mq = 2 TeV with c = 0.1, but the identification would 
be allowed only for Mq = 0.7 — 1.4 TeV with c not larger than, say, 0.05. 
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VI. FINAL CONSIDERATIONS 



In the previous sections, we have discussed the features of different-spin s-channel ex- 
changes in the inclusive diphoton production process ([I]) at the LHC. Specifically, we have 
considered the hypothesis of the spin-2 RS graviton excitation exchange as the source of an 
eventually discovered peak in the diphoton invariant mass, and compared it to the interpre- 
tation of the same peak as a spin-0 scalar exchange exemplified by model J^. The aim has 
been of determining quantitatively the domain in the RS graviton mass Mq and coupling 
constant c = k/Mp\, where the former hypothesis can be identified against the latter (the 
so-called identification reach). 

Clearly, to this purpose, in the situation of equal number of peak events from the two 
models, the information from the distinctive photon angular distributions is needed. The 
relevant cross sections, differential and angular integrated, at next-to-leading order in QCD 
have been introduced for both the RS and the scalar-exchange model. The comparison with 
the leading order calculations shows that the NLO effects are substantial, and non-negligible 
for this analysis. As shown in Sec. IV, in the TeV resonance mass range of interest here and 
at the considered LHC energies, K-factors turn out to be large, of the order of 1.5 or so for 
the RS exchange, and even larger, of the order of 3 or so, for the scalar exchange model. 
Moreover, while K-factors exhibit an angular dependence in the case of the RS model, they 
do not noticeably alter the flat shape of the angular distribution for the pure scalar resonance 
exchange derived at leading order in QCD. 

The angular analysis to discriminate the RS from the scalar model has been based on 
the center-edge asymmetry Aqe, also estimated at NLO in QCD, for both the 14 TeV and 
the 7 TeV LHC. The numerical results for the predicted total number of resonant events 
and the minimum number of events for RS identification, obtained by the simple statistical 
arguments outlined in the previous section, are presented in Figs. [12] and [131, and summarized 
in Table [H 

As regards the 7 TeV LHC, the theoretical results obtained above show that, taking into 
account recent experimental limits on the RS graviton mass, for luminosity around 10 fb^ 1 
there may still be the possibility to identify the RS graviton excitation, in a rather limited 
range of Mq and small c. Higher luminosity, and/or larger LHC energy, would be required 
to extend this region by the Aqe based angular analysis. 
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~0.9 


0.05 
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1.4 


0.1 


4.2 


3.2 


2.0 
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TABLE I: Theoretical discovery (5<r-level) and 95% CL identification limits (NLO) on Mq in TeV. 

Finally, one may observe that the results regarding the RS graviton identification ob- 
tained in the previous section show that the inclusive diphoton process (CQ) at the LHC 
is complementary to the Drell-Yan dilepton production, in the sense that only the scalar 
exchange needs be considered as a competitor, alternative, hypothesis for the source of reso- 
nance events. We have performed the Ace angular analysis of the diphoton production cross 
sections for both hypotheses at NLO in QCD, and estimated numerical results for discovery 
and identification. The Ace method has so far been applied to the dilepton channel at the 
LO in QCD only. Therefore, it should be interesting, for a fully exhaustive comparison of 
the results achievable from the two channels, to extend the AcE-based angular analysis at 
NLO also to the dilepton production process. 
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